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We show how to generahse the zero temperature Lanczos method for calculating dynamical cor- 
relation functions to finite temperatures. The key is the microcanonical ensemble, which allows us 
to replace the involved canonical ensemble with a single appropriately chosen state; in the thermo- 
dynamic limit it provides the same physics as the canonical ensemble but with the evaluation of a 
single expectation value. We can employ the same system sizes as for zero temperature, but whereas 
the statistical fluctuations present in small systems are prohibitive, the spectra of the largest system 
sizes are surprisingly smooth. We investigate, as a test case, the spin conductivity of the spin-1/2 
anisotropic Heisenberg model and in particular we present a comparison of spectra obtained by the 
canonical and microcanonical ensemble methods. 



I. INTRODUCTION 

The study of lattice quantum many body systems by 
the exact diagonalisation technique has proven popular 
at zero temperature (T = 0) where only the ground state 
is required, but it is of less use at finite temperature. The 
reason can be attributed to the different system sizes ap- 
plicable, where for spin-1/2 the ground state can be found 
for up to 7V~30 lattice sites, but the entire spectrum can 
readily be achieved for systems only up to N^IQ. The 
intrinsic difficulties associated with applying the finite- 
size scaling method on such small systems severely limit 
finite temperature applications. At T = the contin- 
ued fraction techniquc^'^ allows accurate calculations of 
dynamical correlations using only the machinery of the 
Lanczos algorithm, but unfortunately this technique has 
not been extended to finite T where mostly full diago- 
nalisation has been employed. As well as direct appli- 
cations of the canonical ensemble, there is also a hybrid 
method which employs the canonical representation of 
dynamical correlation functions but uses a Lanczos basis 
to provide a set of orthogonal states^. This method al- 
lows access to larger systems than are accessible to full 
diagonalisation techniques but to smaller systems than 
the current proposal, which does not need details of all 
the states even in the Lanczos basis. In this article we 
extend the T = formalism to finite temperature by 
applying a microcanonical ensemble approach combined 
with the Lanczos method (MCLM) that provides smooth 
predictions for dynamical correlation functions at least at 
high temperatures. 

The physical advance is to appreciate that in the ther- 
modynamic limit the microcanonical ensemble is equiv- 
alent to the canonical ona^, but for finite systems this 
is much easier to work with. The statistical fluctuations 
engendered by the microcanonical choice are a drawback 
for small systems but become controllable for large sys- 
tems. In practice, as the finite T calculations are much 



smoother, it is more natural to contemplate applying 
finite-size scaling than for their T — Q counterparts. The 
exponentially dense nature of a many particle spectrum 
in the bulk is the property that smoothes our calcula- 
tions, a characteristic that is lost near the ground state 
where the spectrum is sparse. 

Besides the computational interest of this proposal it 
is worth pointing out that, to our knowledge, no studies 
of the fundamental equivalence between the microcanon- 
ical and canonical ensemble for quantum dynamic corre- 
lations exist in the literature. Thus this work is a step in 
numerically exploring this basic postulate of nonequilib- 
rium statistical mechanics; clearly, analytical studies are 
needed to clarify, for instance, the meaning of the mi- 
crocanonical ensemble for a quantum system with dense 
spectrum as an average over a single quantum state (or 
a narrow window of states) and the finite size corrections 
inherent in this ensemble. 



II. THERMODYNAMIC ENSEMBLES 

In this section we will discuss how an arbitrary proba- 
bility distribution can be used, under reasonable assump- 
tions, to represent the canonical ensemble in the thermo- 
dynamic limit. The choice of a (unnormalized) distribu- 
tion p{e) in the thermodynamic limit can be examined 
by considering its Laplace transform. 



e/(^) = de e'^'pie), Kr < 
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where /(r) controls the properties of a distribution des- 
ignated by p(e). As examples, the following choices of 
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/(r) lead to, 

e/(-) = -T^ - p{e) = e-P^ (2) 
the canonical ensemble, 

e/(") = e^" ^ p(e) = 5{\ - e) (3) 
the microcanonical ensemble at energy A, 

the "Gaussian" ensemble at energy A and width a. 
If we examine the partition function Zp, then. 



III. DYNAMICAL CORRELATIONS IN THE 
MICROCANONICAL AND CANONICAL 
ENSEMBLE 

The usually studied quantities of direct physical inter- 
est are the dynamic structure function, 

/+00 
dte*"*(Xq(t)X_q(0)), (8) 
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and dynamic susceptibility, 
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X(q,c.) = I / dte-*([Xq(t),X_q(0)]), (9) 

JO 

where z = lj A- irj, the angle brackets denote a canoni- 
cal ensemble thermal average and the commutator plays 
a central role in the linear response theory (or Kubo) 
formulation of transport. 

The two quantities are related by the fluctuation- 
dissipation relation, 
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where Z[t) is the canonical partition function. The phys- 
ical idea behind the thermodynamic limit is that the par- 
tition function becomes immensely sharp when consid- 
ered as a function of state space and becomes dominated 
by the large number of states with the correct thermody- 
namics; in practice fluctuating quantities can be replaced 
by their thermodynamic average with negligible error. 
Mathematically, an integral such as JSJ may be approxi- 
mated in the asymptotic thermodynamic limit using the 
idea of 'steepest descents' with negligible error. 
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Here fi* is chosen so that. 
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where /3 = l/ksT is the inverse temperature. Note 
that 5(q, lo) satisfies the symmetry relation S'(— q, —u) = 
e"^^ S{(\,uj) while the sum-rule, 
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(11) 



makes it natural to consider the normalised to a unit area 
correlation function, 



(12) 



We have presented the dynamical correlation functions 
in the canonical ensemble and now we will establish their 
form in the microcanonical one. Starting from equation 
© and employing solely the idea that our distribution 
has a restricted energy A we can generate a correlation 
function s(q, w) in the microcanonical ensemble, 



and the average energy at the desired temperature is cru- 
cial. For the particular case of the microcanonical distri- 
bution A = {H) and we need to employ states whose en- 
ergy is the thermodynamic average as one might naively 
guess. Provided that /(r) has only a weak dependence 
on the system parameters, then the two partition func- 
tions are essentially equivalent. It is also clear that pro- 
vided /(r) has the required properties that the 'steepest 
descents' is a good approximation, then any appropri- 
ate ensemble will provide the thermodynamic limit, for 
a particular temperature. It is quite natural to employ 
/(r) = Xt -\-F{t) where A is extrinsic and F{t) is intrin- 
sic, in order to limit towards the microcanonical ensem- 
ble. 
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Here, we have used the relation, 

{OU{H)) ^ {0)U{X) (14) 

{U{H) a function of H) and a decomposition using the 
eigenbasis | m). The expression (|13l) integrates to pro- 
vide, 

s(q,u;) = 27r^(Xq | m){m \ X^^)5{lo + A - £„), (15) 



3 



that can be re-represented as the basic correlation in the 
microcanonical ensemble, 

s(q,a;) = -2 lim 3(Xq [z - + A]"^ X_q). (16) 

Notice that this expression is exact in the zero tempera- 
ture limit where the expectation value is to be taken over 
the ground state wavefunction. 

Now let us imagine that we could find a single eigen- 
state at will, with an energy arbitrarily close (in the ther- 
modynamic limit) to a target energy, A say. It is in prin- 
ciple straightforward then to determine 

s*(q, w) = -2 lim 3(* \ X„[z - H + e,]"^ | *), 

7/1— >0 

(17) 

exactly as before, where _ff | *) = | *) is the known 
eigenstate with e* i— > A. If the microcanonical ensemble 
is equivalent to the canonical ensemble and if a single 
eigenstate is representative of the microcanonical one, 
then provided that A = {H) for the desired temperature, 
we can expect that 

s*(q,w) S'(q,w) (18) 

in the thermodynamic limit. This amounts to the physi- 
cal idea behind our calculations. 

Furthermore, from (|10|l . it follows that 

X"(q,^)-i^|-^s(q,w) (19) 

and from the symmetry of S'(q, lo) in the canonical en- 
semble we can deduce that, 

In^^^p., (20) 

this relation then provides an alternative, cross-checking 
technique for determining the temperature for a partic- 
ular value of A. Although we might like to believe that 
a single eigenstate corresponds to the microcanonical en- 
semble, based on a putative ergodicity assumption for the 
eigenstate, in practice it is not possible to find such an 
eigenstate. So we relax the eigenstate hypothesis and go 
back to a distribution of eigenstates close to the desired 
value A. We simply use the formalism as though we had 
such an eigenstate. 

IV. THE LANCZOS METHOD 

In principle we must construct a particular eigenstate 
with energy A that equals the canonical expectation value 
of the energy, < H >= A, at the desired temperature. In 



practice, we employ the well known Lanczos algorithm 
that is an efficient way of diagonalising large Hamiltoni- 
ans using as variational subspace (truncated basis) the 

set of states, 

{\0),H\0),...,H"^ \0)}, (21) 

where | 0) is a (usually random) initial state and Mi + 1 
the number of Lanczos steps. To obtain an eigenstate 
close to energy A one might expect to use the closest 
eigenstate to A in the truncated basis, but this is to- 
tally incorrect. In practice, only the states at the edge of 
the spectrum converge and the other 'eigenstates' in the 
truncated subspace have the suggested energies but are 
usually far from eigenstates. 

In order to apply the Lanczos method idea, one can 
simply push the energetic region of interest to the edge of 
the spectrum by choosing an appropriate new operator. 
One natural choice is to use, 

K={H-\f (22) 

which is positive definite and pushes the eigenstates with 
energy close to A towards the minimal, zero, eigenvalue 
of K. Another way to understand this technique is to 
consider expanding the ground state of that we will 
call I A), as a probability distribution over eigenstates. 
Choosing A establishes the appropriate mean for this dis- 
tribution but minimising K corresponds to minimising 
the variance of the distribution, and consequently local- 
ising the distribution near A. 

One can perform a Lanczos calculation based upon the 
operator K or, more efficiently, one can evaluate the oper- 
ator K (now a pentadiagonal matrix) in a previously con- 
structed Lanczos basis using H ("L-projection" method). 
Note that, 

{{H {H)f) = {{H Xf) m Xf > 0, (23) 

(the expectation value is over | A)) and so a small variance 
guarantees a narrow distribution of energies around A. 

In any Lanczos calculation the mathematical orthog- 
onality between states becomes lost at some stage as 
numerical errors build up. In practice only the well- 
separated converged states suffer from this disease and 
for us these states, which are at the edge of the spectrum, 
do not gain any significant weight in the correlation func- 
tions and so do not manifest in our results. The states 
at low frequency are all well behaved and maintain their 
orthogonality. 

It is straightforward to implement these ideas numer- 
ically, with a 'double-Lanczos' calculation; the first run 
through a Lanczos procedure of Mi steps is employing 
the operator K starting from a random state and it is 
used to find the state | A) which plays the role of the 
microcanonical distribution. The second run of M2 steps 
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through Lanczos is made using Xq | A) as the initial 
state and then the resulting tridiagonal matrix can be 
diagonalised to form the dynamical correlations directly 
or by employing the continued fractions method which is 
numerically more efhcient but introduces a loss of reso- 
lution. 

All the analysis so far has been subject to several 
caveats; firstly, that the microcanonical ensemble is 
equivalent to the canonical one in the thermodynamic 
limit and in the context of quantum dynamic correla- 
tions. Secondly, that a single eigenstate is equivalent 
to the microcanonical ensemble and thirdly that we can 
find such an eigenstate at will. The first two assump- 
tions, as we have mentioned in the introduction, should 
be the focus of analytical studies as fundamental issues 
of nonequibrium statistical mechanics. 

Regarding the third assumption, it is clearly prob- 
lematic as it is well known that although the Lanczos 
method converges quite easily at the sparse edges of the 
spectrum, in the denser inner regions of the spectrum, 
of interest at finite temperature, it takes the Lanczos 
procedure an exponentially large number of iterations to 
converge. A many-body spectrum has an exponential 
number of states, e.g. for spin-1/2 the # (States)'^ 2^, 
and for a bounded Hamiltonian the eigenstates are com- 
pressed into an energy region that grows only linearly 
with system size. Although the low energy region main- 
tains a sparse density of states, the eigenstates become 
exponentially close together in the area of interest and 
essentially become unattainable. 

At first sight this appears an insurmountable difficulty, 
but in practice this issue allows the technique its success. 
The first Lanczos procedure provides a single quantum 
state I A), that is not an eigenstate, but which when 
decomposed in an eigenstate basis, it is represented by a 
narrow distribution |a„p around A; 

I A) = ^ a„ I n), H \n) = en\n) 

n 

(A I i/ I A) = A, (24) 
gives for the expectation value of an operator O, 

(A I O I A) = ^ I a„ p (n 1 O I n) 

n 

-I- al.anim \0\n). (25) 

n^rn 

This state, used in the evaluation of expectation val- 
ues, acts as a statistical average over an energy window. 
It is important to notice, that by employing a single 
quantum state (not eigenstate) for evaluating an expec- 
tation value (as a substitute for a statistical average over 



a narrow energy window of eigenstates), we assume that 
the appearing off-diagonal terms (second term in ea. (|25(l ') 
cancel each other. This assumption can be justified (and 
numerically verified) by invoking a random phase decom- 
position of the used quantum state. 

From this discussion we can expect two types of fluctu- 
ations in the obtained spectra; first, intrinsic fluctuations 
due to the finite size of the system, present even when 
a single eigenstate is used for the evaluation of the ex- 
pectation value. Second, statistical fiuctuations entering 
by the off-diagonal terms in eq. 1)25(1 due to the use of 
a single pure state that is not an eigenstate; this type of 
fiuctuations can be reduced by averaging over orthogonal 
states I A) (e.g. corresponding to different translational 
symmetry k~ subspaces as we will show below). 



V. CONVERGENCE OF PROJECTION 

In the following we present a test on the rate of con- 
vergence of the projection to a single quantum state with 
energy close to A. Due to the innate complexity of an 
implicit scheme like Lanczos, we develop the theory of a 
simpler technique bricffy to exhibit the ideas. 

A rather simple method of numerically solving for the 
ground state is by an iterative sequence of applications of 
the scaled Hamiltonian. For us this amounts to iterative 
applications of the operator, 



P =1- 



H ^\ 



(26) 



where fi is chosen to be large enough so that /i^ > 
(e„ — A)^ for the full spectrum. Repeated applications of 
this operator exponentially suppresses all states except 
those for which e„ ~ A which remain unaffected. We can 
start out with a set of random states and then for M 
applications of our operator we can build a distribution. 
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{0{P) is the step function) and if we were to perform an 
average over an orthogonal basis, | ^l)), then this would 
converge to. 



Pm{H) ^ 9 [P] P'''. 
Elementary analysis provides: 
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In the limit that M i-^ cxd wc find that, 
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and we converge to a narrow Gaussian probabifity distri- 
bution, 



Pm{H) ^ exp 



(31) 



The width of this distribution is under our control, 



{{H-\f) 



4M + 3 4M + 3 



(32) 



where W is the natural energy scale for the model and we 
see that M needs to scale with the square of the system 
size N to maintain resolution. 

The Lanczos method is clearly much more sophisti- 
cated and provides a much narrower distribution. We 
have examined the distribution obtained in a Lanczos 
calculation and we find that it is well represented by a 
Gaussian distribution with a variance controlled by the 
'eigenvalue' of K attained by the calculation. In prac- 
tice this is about two orders of magnitude better in en- 
ergy than the result obtained from the projection anal- 
ysis fea. (|26|l ). which however it is analytically control- 
lable; indeed, we find that the Lanczos method scal es as 
{K) cx Mi~^ so that the intrinsic resolution, a = \/ {K), 
is inversely proportional to the number of iterations. The 
convergence properties of the three schemes we discussed 
are depicted in Figure ^ for a representative calculation 
of the study that we present in the next section. 

The application of the technique should now be trans- 
parent; employing a single random state, or averaging 
over a sequence of orthogonal random states, one per- 
forms a first Lanczos calculation of Mi steps to find 
the approximate ground state | A) for the operator 
K — {H — Xf. The value of A must be pre-selected 
so that A = {H) for the chosen temperature; several 
techniques are available for reliably determining this en- 
ergy versus temperature relation as the Bethe Ansatz 
(for integrable systems), the finite temperature Lanczos 
(FTLM)^, the Transfer Matrix Renormalization Group 
(TMRG) or Quantum Monte Carlo method. The degree 
of convergence can be measured using the eventual 'eigen- 
value' of K] it plays the role of the variance of the cho- 
sen distribution and its square-root is an intrinsic energy 
resolution a. This scale, cr, can never drop below the dis- 
tance to the nearest eigenvalue. For a usual size system. 
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FIG. 1: Convergence properties of different Lanczos projec- 
tion procedures: (i) dashed line, using eg. 12611 . (ii) dotted line, 
using K = {H — , (iii) continuous line, "L-projection" (see 
text). 



e.g. TV > 16, and temperature, this limit is unattainable 
but a resolution of cr ^ 0.01 {{K) - 0.0001) is readily 
available with a thousand or so Mi iterations. 

Once one has found this state | A), that plays the role 
of the state | *), a second Lanczos projection sequence is 
generated employing the state Xc^ \ A) as the initial state. 
The resolution of the eventual result is controlled by the 
intrinsic dependence on the microcanonical ensemble and 
the degree of convergence measured by a. This can be 
seen from relation (|16|l as the eigenstates over which the 
state I A) is decomposed have a spread in energy a with 
respect to the reference energy A. The resolution also 
depends on the convergence achieved in the 2nd Lanczos 
procedure where the number of iterations M2 denotes the 
finite number of poles which are used to try to represent 
the dynamical correlations. At the sparse edges of the 
spectrum these poles denote the eigenvalues of the sys- 
tem but in the bulk of the spectrum, when grouped into 
bins of a given frequency width, they are fairly uniformly 
spread and offer a further natural energy resolution for 
the calculation. 

More Lanczos steps provide more poles and a finer 
spectral 'grid' for the correlation functions, until the 
graininess of the real system is achieved. We have elected 
to use a few thousand poles in our calculations with very 
little improvement obtained by increasing this number 
as we shall see. The final resolution is self-imposed and 



6 



is the r] of (|16() which we choose to be of order of the 
spectral grid in order to smooth our calculations. 



VI. APPLICATION ON THE SPIN-1/2 
HEISENBERG MODEL 

We are now in a position to test our proposed technique 
and uncover its strengths and weaknesses. We have cho- 
sen to investigate the finite temperature dynamics of the 
prototype spin-1/2 Heisenberg model (equivalent to the 
fermionic — y model). This choice was dictated by 
its central role in low dimensional quantum magnetism; 
an exact solution of the thermodynamics and elementary 
excitations is known using the Bethe Ansatz method^, 
the spin dynamics probed by NMR is of current exper- 
imental interest^ and several numerical and analytical 
studied have been devoted to the study of finite temper- 
ature dynamic correlationsSiiLiSiiiii. The Hamiltonian 
is given by, 
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FIG. 2: Microcanonical calculations for N = 26, A = 2, 
n = 0.02; (a) T = 0, (b) /3 ^ 0. 
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H^Y.^i = jj2isrsr+, + sfs^, + Asfs^,), (33) 



where S"" {a 



x,y,z) are spin-1/2 operators on site I 
and we take J as the unit of energy and frequency {h — 

In particular, we will look at the high temperature spin 
conductivity in the antiferromagnetic regime, J, A > 0, 
for which several studies exist and some exact results are 
known-iS?". To discuss magnetic transport, we first define 
the relevant spin current, , by the continuity equation 
of the corresponding local spin density S[ (provided the 
total component is conserved). 



5^ = ^5f, 



Vjf = 0. 



(34) 



Thus, we obtain for the spin current j^, (that plays the 
role of the operator Xq), 

f = E = JY^^SfS^^, - SfSr+,). (35) 



The real part of the "spin conductivity" (t'(w) (corre- 
sponding to the charge conductivity of the fermionic 
model) includes two parts, the Drude weight D and the 
regular part (Treg(w)24ifi, 



lj'{u!) = 2TrDS{uj) + Gregi^)- 

The regular contribution is given by. 



(36) 



1 — e~^'^ n v-^ 

CTres(w) = — y ^ Pn\ < n\f\m > ^ (5 (w - LJ„i„ ) 



N ^ 



(37) 



where p„ are the Boltzmann weights and = Cm — 

To compare the presented data on the conductivity we 
normalize them using the well known optical sum-rule 
that in the /3 ^ limit takes the form, 



-t-oo 



(38) 



The normalized conductivity, cr(i^), in this high temper- 
ature limit is given by. 



cr(w) 



if') 



(39) 



that can be calculated using our microcanonical ensemble 
procedure by. 
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(40) 



In principle this expression includes also the contribution 
from the zero frequency Drude weight (5— function, but 
in practice as the second Lanczos procedure cannot fully 
converge, the Drude peak appears as a low frequency 
contribution. As we will discuss below, sorting out this 
low frequency part, in general allows us to reliably extract 
the Drude weight value. 

In general, we can employ the translational symmetry 
of the Hamiltonian and study spectra in a given pseu- 
domomentum k— subspace or average the results over 
different k— subspaces; in the following we typically em- 
ploy Ml = 1000 and M2 = 4000 Lanczos iterations at 
/3 — ^ unless otherwise stated. In FigureElwe compare a 
zero temperature^ with an infinite temperature {(3 — > 0) 
calculation for a fairly large system in the fc = subspace. 
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FIG. 3: Microcanonical versus Canonical calculations; (a) 
iV = 20, A = 0.5, Tj = 0.01, (b) iV = 18, A = 1, 77 = 0.01. 



The zero temperature calculation finds a few poles with 
exact weights whereas the infinite temperature calcula- 
tion provides a much smoother result. 

There is clear structure in the infinite temperature re- 
sult but also apparently some noise. To interpret this 
result we must consider the issue of the veracity of the 
microcanonical ensemble for such small systems namely 
the extent to which the microcanonical ensemble is equiv- 
alent to the canonical one. 

In Figure |31 we present a comparison, extremely en- 
couraging, of some microcanonical calculations with the 
analogous canonical ones. There is 'noise' in all calcu- 
lations, the origin and magnitude of which we will now 
discuss. The canonical calculations are essentially a di- 
rect evaluation of expression 137|l . where we applied a 
"binning" procedure on the (5— function weights over an 
energy scale of about 0.01. The number of contribut- 
ing matrix elements are of the order of the dimension 
V of the Hilbert space squared, 2?^, e.g. 10^ — 10^ S— 
functions, with no continuity in the weights. The results 
are not smooth and the resulting intrinsic fiuctuations 



are heavily smoothed by our binning procedure. In the 
microcanonical calculations we employ our scheme, fur- 
ther averaging over translational symmetry fc— subspaces. 
Now, only of 0{'D) 5— functions are essentially contribut- 
ing, multiplied by the number of states involved in the 
decomposition of the state | A) (a few thousand depend- 
ing on the convergence) and the number of fc— subspaces. 
We could average over initial random states, but we find 
that this has only a small smoothing effect, because the 
underlying poles are at the same energies. Notice that 
the observed fluctuations are not associated with any of 
our different resolution processes which are much smaller 
than the observed scale of fluctuations; they are due to 
the finite size of our system and thus to the effective 
smaller number of matrix elements contributing to the 
construction of the spectra. This seemingly new problem 
associated with our technique turns out to be dominant 
for small system sizes; very soon however it becomes neg- 
ligible as larger systems are achieved, specially consider- 
ing that the dimension of Hilbert space grows exponen- 
tially fast with the system size N. 

In order to assess these fiuctuations and simultane- 
ously the role of our smoothing parameter 77, we per- 
formed some basic calculations involving only a single 
fc— subspace state | A). In Figure^lwe offer a comparison 
of calculations involving just the poles evaluated using 
the 2nd Lanczos procedure eigenstates against smoothed 
versions of the same data but employing the continued 
fraction technique. 

The fluctuations clearly decay with system size with 
the flnal system being surprisingly smooth. The limi- 
tations of the smoothing process are clear, the sharper 
features are slightly washed out although the ease of as- 
sessing the data makes such a smoothing advisable. The 
weights for these microcanonical calculations are truly 
quite continuous in comparison to the intrinsic properties 
of the canonical calculation which is necessarily ragged. 
Obviously for our largest calculations we are nowhere 
near converged to the true spectrum which is a possible 
explanation for the observed continuity. 

We can now fairly safely conclude that our technique is 
a viable way to calculate dynamical correlation functions 
at high temperature for the same systems accessible by 
the Lanczos method at T = 0. By its very nature, the 
finite T correlations are much smoother and more regular 
to interpret. Our technique introduces new statistical 
fiuctuations which make small system sizes ragged but 
appear to leave large system sizes essentially unaffected. 

Although we can now investigate finite temperature 
dynamic correlations using the Lanczos method, we are 
still restricted to iV ^ 30 for a spin-1/2 system. The key 
to making useful physical deductions is the procedure of 
finite-size scaling, the attempt to deduce the properties 
of the infinite size system using assumed properties of 
the size, N, dependence. This method has been widely 
and succesfully applied in the evaluation of ground state 
energies or gap values using data provided by the exact 
diagonalization, Lanczos or Density Matrix Renormal- 
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ization Group technique. But to extract information on 
finite temperature dynamic correlations one would need 
to know the form of the curves before fitting and scal- 
ing could take place mathematically. As it is clear from 
Figure [SI this might be a challenging task considering the 
statistical fluctuations inherent in the spectr however, 
from ongoing studies on other systems using this method, 
we find that the behavior of the spectra might greatly de- 
pend on the model Hamiltonian and correlations under 
study (e.g. it is far more structurless for energy current 
dynamic correlations). Note that the high frequency be- 
havior is generally rather weakly size dependent while the 
low frequency one is the most subtle to determine. The 
last however is the most physically interesting as it de- 
termines, for instance, the diffusive or ballistic behavior 
of the conductivity. 

The basic properties of the /3 — > current-current cor- 
relations are now available and so we provide in Figure 
a few examples of the frequency dependence of the con- 
ductivity at /3 — > as a function of A. 

Although we have devoted most of our effort to infi- 



FIG. 6: Microcanonical ensemble evaluation of the normal- 
ized conductivity a{uj) for /3 ^ 0, = 28; (a) A = 0.5, 1.0, 
(b) A = 2.0,4.0' 

nite temperature (/9 — > 0), our technique is valid at es- 
sentially any temperature (provided that we remain at 
a dense region of the spectrum). Analysing the pure 
Heisenberg model, we look at a couple of finite temper- 
ature fc— averaged calculations in Figure The temper- 
ature has been deduced from a least-squares fit of the 
quantity, 

In a + Pmicro^ (41) 

to a linear Ansatz, and although the statistical fluctu- 
ations are compounded, an almost vanishing intercept 
and a clear slope indicate the feasibility of the strat- 
egy. The obtained Pmicro values compare favorably with 
those corresponding to the canonical ensemble in the 
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FIG. 7: Finite temperature calculations for A'^ = 24, A=l, 
T] = 0.01; (a) A = —3, (b) Temperature fit fimicro — 0.14, 

(c) s{lo), X = —6, (d) Temperature fit jSmicro — 0.28. 



thermodynamic limit, evaluated using A =< H >p; for 

A = -3, Prmcro ~ 0.14 VS. Pcanomcal ^ 0.15, for A -6, 
f^micro ^ 0.28 VS. Pcanonical ^ 0.3. 

Although we have compared numerical evaluation of 
dynamic correlations obtained by a canonical and micro- 
canonical method, we have yet to compare with an exact 
solution. Recently even non-zero temperature dynamical 
correlations have become partially accessible, with a cal- 
culation of the Drude weight for the < A < 1 Heisen- 
berg model at finite temperaturei^. In particular, the 
Drude weight in the /3 ^ limit is given analyticalljiii 

by, 




^ ^ l (V^-0.5sin(2V^)) 



cos(7r/z/). (42) 



FIG. 8: A comparison of three 'microcanonical' distributions 
(K) = 0.002(Mi = 500),0.0005(Afi = 1000), 0.0012(A/i = 
2000), for A = 26 and A = 0.5; inset, low frequency range. 



The Drude weight, strictly speaking, is defined as the 
weight of a zero frequency (5— function, ea. (|36|l : it is a 
particularity of the Heisenberg model that it appears as 
a narrow peak at low frequencies, of the order of the 
inverse lattice sizei&, in contrast to the fermionic "t-V" 
version where it is accounted for only by the diagonal 
energy elements (w = 0). 

In extracting the Drude weight by the above described 
procedure we must take into account the problem caused 
by the intrinsic resolution of our calculations, which is 
of order a = y/{K). Although our chosen resolution of 
a ~ 0.01 is almost invisible for the smooth background, 
for the Drude weight the resolution is essentially limited 
by that of our 'microcanonical' distribution, viz a. An 
example of these ideas is provided in Figurc|Hl from which 
it is clear that the Drude peak is the only contribution 
for which the change in resolution is relevant. These 
calculations involve a single state and are much improved 
by fc— averaging, also the energy window is so small that 
the individual poles in the 2nd Lanczos procedure are 
visible and have been smoothed out with an 7y=0.005 
which adds to the observed resolution. In the inset, the 
scale of the conductivity clearly signals a low frequency 
peak (notice the difference in scale between Figure |H1 and 
its inset); still in order to extract the Drude weight from 
the smooth background, we must integrate the peak up 
to at least as far as it is resolved and that necessitates 
the inclusion of some of the background. We have elected 
to err on the side of inclusion and tend to integrate past 
where the Drude peak appears to become small. 

In Figure El we offer a comparison of the analytical 
and numerically extracted Drude weights in the /3 — » 
limit. The quantitative agreement is reasonably satisfac- 
tory, becoming rather poor near A ~ 1 because of our 
technique for extracting the Drude weight; due to the 
finite resolution of our calculation we need to sample a 
finite width around w = 0. For the case A = 1 there is 
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FIG. 9: Comparison of /3 ^ Drude weight, D / P; numerical 
evaluation (points) vs analytical expression eq. 1421 (contin- 
uous line). 



no Drude weight but there does appear to be a power- 
law like divergence which we pick up in our finite window 
leading to the observed corrupted behaviour. 



VII. DISCUSSION 

Our investigation appears to validate the use of the 
Lanczos algorithm to analyse finite temperature dynam- 
ical properties of strongly correlated systems; the crucial 
step is to employ the microcanonical ensemble, which 
essentially allows the thermodynamic average to be re- 
placed by an elementary expectation value. All the sim- 
plicity of the zero temperature formalism can then be 
taken over to the finite temperature calculation. The 
comparison of canonical with microcanonical procedures 
indicates that the thermodynamic limit is reached with 
quite modest system sizes and consequently there ap- 
pears to be little systematic error coming from our choice 
of ensemble. There are intrinsic statistical fluctuations 
in our calculations but these are severely curtailed by 
increasing the system size and are an implicit difficulty 
with canonical calculations too. We believe that we can 
calculate the high enough temperature dynamical corre- 
lations for a finite system with an excellent tolerance. 

The statistical fluctuations in our results require to be 
controlled if an error analysis is to be contemplated. Al- 



though we have not got analytical control, we do have 
experience at various approaches to reducing the statis- 
tical fluctuations. The crucial point is that, when taking 
a statistical average, one should use "orthogonal" states 
(I A}'s decomposed into different sets of eigenstates | n)). 
Averaging over random starting vectors in the same sub- 
space is not very effective, even if they are originally or- 
thogonal, because the resulting distribution involves the 
same states and consequently an overlap. Performing a 
k— average, over translational symmetry subspaces, is an 
excellent procedure, since the states are automatically or- 
thogonal and intellectually one is reverting back towards 
the real physical statistical average. Another possibility 
is to use several of the eigenstates of the first Lanczos pro- 
cedure; although the orthogonality is guaranteed, there 
is an induced loss in resolution due to the larger cr's of 
the higher Lanczos states. A final possibility is to em- 
ploy the parameter A, where the average over different 
A's must be limited within a window that corresponds 
to the energy fluctuations at the studied temperature in 
the given size system. Providing that the A's are further 
apart than the chosen a, the orthogonality is essentially 
guaranteed. 

Although we believe we have access to the temperature 
behaviour of finite-size systems, this does not give imme- 
diate access to the dynamics in the thermodynamic limit 
because finite-size scaling must be performed; Figure 6 
exhibits clear peaks of unknown form, plausible 'cusps' 
and regions where the correlations vanish. Unless we can 
guess or deduce the form of these structures, finite-size 
scaling appears problematic. We should note however 
from our experience, that not all models and dynamic 
correlations exhibit so involved spectra; in forthcoming 
works we will present analysis of charge/spin/energy cur- 
rent correlations for other (non-) integrable systems of 
current interest (higher spin, ladder models) where the 
obtained spectra are far more structrurless. Finally, be- 
sides the finite frequency behavior, our method allows the 
reliable study of scalar quantities as the Drude weight. 
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